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'k>( ! Abstract 

j^ ■ In SUGRA inspired supersymmetric models with universal boundary conditions 

for the soft masses, the scalar cross section dscaiar for the elastic neutralino-nucleon 
scattering is in general several orders of magnitude below the sensitivity of current 
experiments. For large tan/? and low Mi/2,mQ values, the theoretically predicted 
o'scaiar Can approach the sensitivity of these experiments (ss 10~^ pb) being at 
the same time in agreement with recent cosmological data, which impose severe 
restrictions on the CDM relic density, and with accelerator experiments which put 
lower bounds on sparticle and Higgs boson masses. Further improvement of the 
sensitivity of DAMA and CDMS experiments will probe the large tan/3 region of 
the parameter space in the vicinity of the boundaries of the parameter space allowed 
by chargino and Higgs searches. 



As time goes by there is indisputable, accumulative experimental (observational) evi- 
dence that we live in a flat and accelerating universe. Indeed, the most recent confirmation 
of such a dramatic statement have been the new results from the BOOMERANG [jl| and 
MAXIMA PI experiments, where for the first time, it has been possible to measure several 
crucial cosmological parameters simultaneously. The position of the first acoustic peak 
of the angular power spectrum strongly suggests that Qq = 1 (i.e. k = 0, a fiat uni- 
verse), while its shape is consistent with the predicted inflationary density perturbations. 
It is rather heartwarming to see the way that the data favour an almost fiat Universe 
Qq = r2jv/ + ^A = 1 ± 0.2 Pi where Qm is contribution of the matter density and Q/^ 
is that of the vacuum energy. At the same time one determines that Qm = 0.4 ±0.1 
which in conjuction with analysis from high-z SNa data result to Q/^ = | Qm + 1=^1 0- 
With Qm = 0.4 ± 0.1 this implies that Q/^ = 0.85 ± 0.2 [|]. Taking into account the 
fact that the baryonic contribution to the matter density is small, Qb = 0.05 ± 0.005, 
the values for matter energy density Qm result to a Cold Dark Matter (CDM) density 
i^cDM — 0.35 ±0.1, which combined with more recent measurements of the scaled Hubble 
parameter ho = 0.65±0.05, result to small CDM rehc densities ^cdm ^o^ — 0.15±0.07 p. 
Combined analysis from MAXIMA, BOOMERANG and COBE/DMR data entails to even 
tighter limits 0, though we are not using such stringent limits trying not to extremely 
constrain the parameter space of the MSSM. 

While the cosmological picture seems to be brighten up at least at the level of a 
phenomenological understanding, we are facing some squeezing at the other end, i.e. at 
the fundamental particle physics front. The question of "who ordered this?" for the 
cosmological constant is so obvious that needs no further elaboration. On the other hand 
the "observed" value of Qm cries out for an non-baryonic component, rather difficult to 
be provided by the nearly massless neutrinos, if we take at face value the present values 
of neutrino masses. 

In i?-parity conserving supersymmetric theories the lightest supersymmetric particle 
(LSP or x) may be a neutralino, which is a good candidate to play the role of dark matter 
(DM) particle. While it is very encouraging that the neutralino parameters (masses, cou- 
plings) are in the right range to provide easily the right amount of DM, e.g. the "correct" 
contribution to Qm, neutralinos still have escaped any direct or indirect detection so far. 
The direct way to detect the neutralino is via the neutralino-nucleus scattering. This 
scattering contains spin-dependent as well as spin-independent (scalar) parts. For the 
heavy nuclei detectors the spin-independent part dominates, because there neutron and 
proton scattering amplitudes are approximately equal. The present sensitivity of experi- 



ments DAMA [^ and CDMS |^ measuring the spin- independent neutralino-proton cross 
section (Tscaian IS approximately 10^^ — 10^® ph. The CDMS collaboration has reported 
negative results from their experiments 0, setting an upper limit on the spin-independent 
o^ scalar that cxcludes also almost all the range suggested by the DAMA collaboration Pj. 
Thus in this letter we focus our interest mainly in the CDMS limits. 

Since neutralinos seem to play a very important link between the very-very small 
(SUSY theories) and the very-very large (DM) it is imperative at every stage to have 
an integrated picture of what is allowed and what it isn't. During the last year we have 
embarked in such an effort P], having used cosmological data and electroweak (EW) 
precision data (mainly from LEP) to constrain, in generic SUSY theories, the parameter 
space. In this letter we continue our programme [§ by injecting the new cosmological 
data and new SUSY bounds into our analysis but take also into account the most recent 
CDMS direct-search exclusion limits. 

As we shall see, in general, one has to improve considerably the sensitivity of the direct- 
search limits in order to start excluding sections of the otherwise allowed parameter space. 
DM detectors have a current sensitivity ~ 10^^ — 10^® ph and would be interesting to find 
regions of the MSSM parameter space where bounds imposed by CDMS experiments are 
nearly saturated. In most of the parameter space the scalar LSP nucleon cross section 



X 



r^^ 



(^scalar IS Small getting values as small as ~ 10^ ph or even less for LSP masses m 
200 GeV. This indicates that DAMA and CDMS data at the current stage are not that 
useful to probing MSSM with universal boundary conditions (UBC) for the soft masses 
at the Unification scale. Relaxing the UBC can result to enhancement of the (Tscaiar as 
has been studied in Refs PHTT . 



Our purpose in this paper is to keep on UBC for all soft masses and focus mainly on 
points that are compatible with neutralino relic densities within the cosmologically allowed 
domain ^^^l ~ ^-^^ '^ 0.07, which yield cross sections agcaiar close to the boundary of 
the area excluded by CDMS experiments. For low tan/3 < 20 this is unlikely to occur 
but for higher values cross sections increase with tan/3 due to the fact that Higgs bosons 



become lighter |jTOl. It is interesting to point out that large tan/5 offers an alternative to 
neutralino-stau coannihilation mechanism to decrease neutralino relic densities, in cases 
the LSP is a high purity bino. Its small Higgsino component is very important since it 
allows it to couple with the pseudoscalar Higgs boson A with a strength that becomes 
sizable for large tan/? despite the smallness of the LSP's Higgsino composition. This in 
combination with the fact that A mass decreases with large tan/3, making pseudoscalar 
exchanges less suppressive, can lead to sizable x X ~^ ^ ^ and x X ~^ 't ^ cross sections 



and hence low relic densities. This has been emphasized in Ref. 0. Although the role 
of the pseudoscalar A is very significant for this mechanism, A itself plays no important 
role in a scalar except that its mass is simultaneously decreased, as tan/3 grows, with that 
of the heaviest of the CP-even Higgses which does mediate the neutralino-nucleon elastic 
cross section. Therefore it would be interesting to perform a refined scan of the parameter 
space, for large tan/?, to search for points which lead to both small relic densities and 
values of ascaiar in the vicinity of the boundary excluded by CDMS data. 

In our analysis we assume UBC for the soft masses and focus mainly to large tan/3 
although results for low tan (3 will be also reported. In our evaluation of the scalar cross 
section a scalar we are using the formulae of Ref. ||T2[. For the hadronic matrix elements 
we use the values of the Ref. [|r^ . 

Regarding the mass of the pseudoscalar boson which plays an important role in de- 
creasing the LSP relic density a subtlety arises related to the theoretical determination 
of its mass which we should now discuss. It is well known that radiative corrections to 
its mass usually calculated through the effective potential approach are not stable with 
changing the scale at which they are calculated. Empirically one calculates its mass at an 
average scale Qt for which the radiative corrections to its mass, due to the third generation 
sfemions, are small and hence can be neglected. In this scheme the pseudoscalar mass 
squared is given by m\ = —2 m^lQ t)/ sin 2 f3{Qt), where r/ig is the Higgs mixing param- 
eter. Although this is in principle correct contributions of charginos and neutralinos are 
not small at this scale especially when M1/2 is large. In some case this may produce an 
error in the determination of its mass as large as 25%. Excluding the chargino/neutralino 
contribution is legitimate provided the relevant scale is not taken to be Qt, but rather the 
average chargino/neutralino mass Qmos- At this scale their contributions can be safely 
neglected. This scale however may differ substantially from Qt-, when M1/2 is large, and 
third generation contributions are not small if they are calculated at Qmos- In our ap- 
proach for the determination of the pseudoscalar's mass we have duly taken into account 
all contributions including the chargino/neutralino corrections, as well as the small gauge 
and Higgs boson contributions, and have observed that they contribute significantly to 
the stabilization of the pseudoscalar Higgs boson mass with changing the scale Q. If these 
contributions were neglected stabilization would be spoiled especially in the large M1/2 
region where the running of the parameter m\{Q) due to gauginos becomes important 
and their corrections should be included to render a Q independent result. The situation 
would be even more dramatic if in addition to having large M1/2 we are in the large tan (3 
regime where such deviations from stability are enhanced as being proportional to tan/5. 



It is also well known that the mass calculated through the effective potential differs from 
the pole mass by n(0) — Il{m\), where n(p^) denotes the corrections to the pseudoscalar 
propagator. For a more reliable estimation of the pseudoscalar's mass we have also in- 
cluded the leading logarithmic parts of 11(0) — Il{m\) to the effective potential mass 
leaving aside unimportant, at least for our analysis, small corrections. In this way we 
have approximated satisfactorily the pole mass and avoid the complexities of calculating 
one-loop integrals, which are usually expressed by Passarino-Veltman functions, as would 
be demanded if we were to determine the pseudoscalar mass through the location of the 
pole of the propagator [|l^ . 



For consistency, both for the calculation of the pseudoscalar mass and for determin- 
ing the Higgsino and Higgs mixing parameters fi, m\ through the one-loop minimization 
conditions, we take into account the contribution of all sectors to the derivatives of the 
effective potential not just those of the third generation. These are given in Ref. ||14 |. 



Since however they are given in the Landau gauge, while we are using two-loop RGE's 
in the 't Hooft's gauge, we remedy this situation by using the gauge and Higgs boson 
contributions to the minimization conditions as these follow directly from the tadpole 
graphs calculated in the 't Hooft's gauge, as can be found in Ref. |T^. More details can 



be found in Ref. 13 



Regarding the neutralino relic density we calculate it in the way prescribed in Ref. 
IQ. We solve numerically the Boltzmann transport equation ignoring at the first stage 
coannihilation effects |jl6||, and in particular f — x and in general / — x coannihilation 



effects |17|. The latter are important for a bino like LSP when its mass is close to that 
of the lightest of the f. The importance of these in reducing the relic density has been 
stressed in Refs [|T7|. One can include the contribution of the coannihilation processes, 
whenever they are of relevance, following the empirical rule 

n^ = R{AM)nl. (1) 

The reduction factor R{AM) depends on AM = {rrifj^ — m^/m^ and the function 
R[AM) smoothly interpolates between k, 0.1 and 1.0 for values of AM in the range 
0.00 — 0.25 (see Ref. |^). The equation above is a handy device and reproduces the 
results cited in Ref. [T^. fi° appearing on the r.h.s. of the equation above is the relic 



density calculated ignoring coannihilation channels. All effects of coannihilation processes 
are effectively included within the function R. 

As a preview of our results, we have found that the largest a scalar i approaching values 
10^^ — 10^^ pb, are obtained in regions of parameter space for which tan/3 is large. The 
dominant contribution to this regime is the Higgs boson exchange. For given inputs 



rriQ, Ml/2, ^0 and the sign of /i, Higgs masses decrease as tan/? increases. The maximum 
value of this angle is determined by the theoretical requirement that we are at the correct 
electroweak minimum and Higgs masses squared are positive and bounded from below by 



the recent experimental limits [jT8| . Hence the contribution of Higgs bosons to neutralino- 
quark elastic cross section becomes more important in the large tan/3 regime. Such a 
decrease in the mass is not sufficient by itself to increase (Jscaiar to levels approaching the 
sensitivity of ongoing experiments. The major role in this increase plays the coupling 
of the CP-even heavy Higgs whose coupling to rf-quark is proportional to ^^. This is 
proportional to tan/5, when the latter becomes large, since cos/3 ~ l/tan/5 and the Higgs 
mixing angle a becomes small and negative in this region resulting to cos a ~ 1. The 
situation alters for the light CP-even Higgs boson whose coupling is ^^. In this case 
sin a behaves as l/tan/5 and unlike the heavy Higgs case its coupling does not grow with 
increasing tan/3 but stays constant of order unity. Therefore despite the fact that the 
heavy CP-even Higgs is heavier than its light CP-even counterpart, its contribution may 
be much larger in the large tan/3 region, due to its enhanced coupling to (i-quark |]I9| . 

From the above argument we conclude that points of the parameter space that are 
likely to yield the highest possible elastic cross sections are those for which tan (3 is large 
and the heavy CP-even state receives its minimum value allowed by experiments and 
other theoretical constraints. The contributions of sfermion exchanges to the amplitude 
of the elastic cross sections are less important and will not be discussed. In the constrained 
superymmetric scenario the mass of the CP-even heavy state mn, has the tendency to 
increase with increasing the effective supersymmetry scale and therefore low rriH values 
are obtained for low mo, M1/2 values. The mass rriH is bounded by the mass of the light 
Higgs boson rrih since 



2 2 I 

rriH = m^^ + 



- {{m\ + m^zY + ^^) "" ''^A "^1 (cos 2/3) + -cos 2/3 {m\ — m|) 



ni/2 

• (2) 



In this expression rriA refers to the radiatively corrected pseudoscalar Higgs boson mass 
and e are the leading stop corrections to the CP-even Higgs masses. Allowing for addi- 
tional non-leading contributions, or for contributions from other sectors, the expression 
above is modified but for the sake of clarity in presenting our arguments we will use the 
simplified relation as presented above. In our numerical analysis we have duly taken into 
account all stop corrections to the CP-even mass matrix as well as those of sbottoms and 
staus. 

In the large tan/? regime cos 2/? ^ —1 and the relation above is simplified to 

ml = ml + -{ml-m%-e). (3) 



Obviously the lowest rriH values are obtained in the region where rriA is light and rrih 
is close to its lower experimental bound. The value of e can also decrease the mass 
rriH, however this depends logarithmically on the stop masses and unlike rriA does not 
vary much with changing the SUSY breaking parameters mo, M1/2 and can be considered 
practically stable. Since Higgs masses increase with increasing the parameters Mi/2,fnQ 
and get lowered with increasing tan/5, the region of interest the most likely to yield 
scalar cross sections nearly saturating the sensitivity of current CDMS experiments is 
the region of large tan/5 and low Mi/2,mo close to the boundaries allowed by Higgs 
boson and chargino searches. It becomes evident that increase of the sensitivity of CDMS 
experiments will probe this region and may exclude points that would be otherwise allowed 
by accelerator experiments. 

It is worth pointing out that in the large tan/? region neutralino relic densities de- 
crease as we have already emphasized due to both the decrease of the pseudoscalar mass, 
whose exchange in xx —>■ bb , rf processes is less suppressive, and the increase of the 
XX^ coupling. The smallness of the LSP's Higgsino component is compensated by the 
largeness of tan/5 yielding neutralino annihilation cross sections compatible with the re- 
cent astrophysical data when M1/2 gets values < 200 GeV [If. Hence there are regions 
in which we can obtain both low relic densities and high agcaiar and these regions can be 
possibly probed by the next round CDMS experiments. 

For our numerical analysis we have scanned regions of parameter space with values 
tan/5 = 1.8 - 50, M1/2 < 1350 GeV, m^ < 1 TeV, \Ao\ < 500 GeV. The dependence 
on Aq is rather mild and for this reason we cut points with values \Aq\ > 500 GeV. In 
Figure |l] we display the behavior of cTgcaiar for both positive and negative values of the 
parameter /i for given niQ = 200 GeV, Aq = GeV and values of M1/2 = 200, 400 
and 600 GeV respectively. One observes that for M1/2 = 200 GeV, which is close to 
the minimum value allowed by chargino searches ( M1/2 )\min ~ 170 GeV, the scalar 
cross section increases with increasing tan/3. Especially in the fi < case this increase is 
more steep and (Jgcaiar can reach values slightly above 10~® pb. In the other cases shown, 
corresponding to M1/2 = 400, 600 GeV, the situation is quite different and Cscaiar gets 
smaller by at least an order of magnitude, due to the heaviness of all sectors involved, 
with the exception of the light Higgs whose contribution in this region is however small. 
The abrupt stop of all lines towards their right endings occurs since for higher values of 
tan j3 we enter regions which are theoretically excluded. Towards the left endings of these 
lines the light Higgs mass approaches its lowest experimental bound ^ 110 GeV. For 
the negative /i case a steep deep is observed for M1/2 = 400,600 GeV and tan/5 = 10. 



This is accidental and due to the cancellation of sfermion and Higgs contributions to the 



amplitudes of the elastic x^nucleon cross section [|ll . 

In Figure |] we plot the scalar cross section as function of the LSP mass m^. Each 
point, struck by either a plus or a cross, has been picked from a sample of 5000 random 
points in the region of parameter space mentioned previously. We simultaneously calcu- 
late the neutralino relic density and we denote points which are cosmologically allowed, 
^x^Q ~ 0.15 ± 0.07, by a plus and those which are not allowed by a cross symbol. 
In this random sample the relic density is calculated without taking into account the 
coannihilation processes. We shall return to this point later. From this figure it is seen 
that irrespectively of the value of the relic density, points that have large cross sections 
^ 10~^ pb correspond to LSP masses m^ < 80 GeV, or equivalently values of M1/2 less 
than about 200 GeV, that is in the region near the edge allowed by recent chargino 
searches. Moreover these points are characterized by large tan/3 (> 30). Some remarks 
are in order. The first concerns the blank stripe occurring at values of my. near 175 GeV. 
This indicates the presence of a top threshold in xx ~^ tt annihilation cross section. This 
appears because in order to speed up calculations we have avoided calculating the relic 
density at singular points, near thresholds or poles, where non-relativistic expressions usu- 
ally employed break down, and more refines techniques must be used to get the correct 
result. Besides this we observe the development of other stripes too in the vicinity of 
~ 100 GeV which correspond to poles of Higgses and other sparticles involved. 

The void region in the negative /i case for values 100 < m^ < 280 GeV, is due to 
the fact that for given value of tan/?, mo and Aq, an upper bound is set on M1/2, or 
equivalently rrij^, beyond which we run into situations where EW symmetry is not broken 
to the correct vacuum. The larger the value of tan f3 is the lower the value of the upper 
bound of Ml/2 is. Hence points of the sample that yield large cross section for low M1/2 
values, for values higher than the upper bound set on M1/2 are completely absent from 
the figure since are automatically excluded. The dispersed points appear for large values 
of Ml/2 in the yU > case are related with the behavior of the (Jscaiar as a function of 
tan/3 for large values of M1/2. Actually for such large values of M1/2 the light Higgs 
boson exchange dominates the XQ ^ XQ process and therefore the (J scalar- However this 
contribution is proportional to the xxh coupling, which varies very much with the tan/3 
for large values of M1/2, resulting to the dispersion shown in the figure. 

Figure |^ follows from Figure |] if we discard points that do not fall within the cos- 
mologically allowed region f^^ ^0 ~ 0.15 ± 0.07. We observe that only a few points 
survive for both negative and positive /i cases when we enforce this strict cosmological 



bound on the relic density. Since we have so far neglected coannihilation processes in our 
analysis, and especially f — x which are the important ones when x is mostly a bino, 
this sample is expected to be enriched when these processes are taken into account. This 
is done in Figure ^ where it is shown that the sample of cosmologically allowed points 
is indeed enriched but not in the region which yields the largest possible cross sections. 
The relic density for points that coannihilation processes are of relevance are calculated 
as explained in the general discussion earlier in this paper. 

In Figures ^ we display a scalar as function of tan {3 for the random points disused in 
the previous figures. It is seen that points in the large tan/5 region approach 10~^ p6, in 
the case /i > 0, and can slightly exceed this value for yU < 0. These are characterized by 
large tan/3 and low mo,Mi/2 values and are in agreement with accelerator experiments 
on sparticle and Higgs searches. Some of these points can survive the strict relic density 
bounds indicating that an improvement of the sensitivity of CDMS experiments can ex- 
plore the parameter space of MSSM with large tan/5 and low mo, Mi/2- Large tan/5 for 



yU > are also compatible with 6 — > s + 7 data from CLEO pO| , and hence next run 
CDMS experiments may be capable of imposing bounds relevant to regions of parameter 
space which accelerator experiments have not probed as yet. 

Concluding we have seen that in order to have as large agcaiar as possible and close 
to CDMS experiment sensitivity, that is 10"^ — 10^^, M1/2 must be as small as allowed 
by chargino searches and tan/5 as large as possible for the Higgs states to be as light 
as allowed by theoretical constraints and experimental searches. This happens both for 
/i > and /i < 0. In addition in the /i > case we can also obtain large ascaiar for large 
Ml/2 ~ 1 TeV, again for large tan/5 > 30, but these points tend to give unacceptably 
large relic densities, fi^ h"^ > 1. On the other hand for /i < and large tan/5, we can get 
large dscaiar for moderate values of M1/2 ~ 450 GeV, which are cosmologically acceptable. 
In all cases considered the necessary condition in order to have ascaiar ~ 10^^ — 10~^ is 
that tan/5 > 35, as can be clearly seen from Figure ^ It is perhaps important to note 
that although the /x < case yields, in general, larger cross sections than the /i > case. 



the former is disfavoured in view of 6 ^ s + 7 experimental data ||20 |. 
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Note added 



After submitting this paper for publication we became aware of the paper |hep-ph/0012377 
by Bottino, Fornego and Scopel in which values for the scalar neutralino-nucleon cross 
section (agcaiar) are obtained which can explain the DAMA data. The region of the 
parameter space in this paper in the case of the mSUGRA model is extended to values of 
mo up to 3 TeV. In our paper we have taken values for the tuq parameter not exceeding 
1 TeV. This range is the same with that explored in the paper |hep-ph/0010203| by 
Bottino, Donato, Fornego and Scopel. The conclusions reached in that paper are in 
agreement with those presented in our work (see figure la of the aforementioned article). 
Increasing the upper limit of rriQ to include values 1 TeV, we enter regions in which 
acceptable neutralino relic densities are obtained with higher values for a scalar^ falling 
within the DAMA experimental sensitivity, at the cost of having a very massive sfermion 
spectrum. The region itiq > 1 TeV with M1/2 < 600 GeV characterizes the "Focus Point 
Supersymmetry" ( see the references J. Feng, K. Matchev and F. Wilczek, Phys. Lett. 
B482 (2000) 388; J. Feng and K. Matchev, Phys. Rev. D63 (2001) 095003). Nevertheless it 
should not escape our attention the fact that the "focus point" region is now dissalowed by 
the BNL E821 experiment (H. N. Brown et ai, Muon g — 2 collaboration, [hep-ex/0102017| ) 
and other constraints as emphasized recently in Ref. Piep-ph/0 102331] on g — 2 (J. Ellis, 
D. V. Nanopoulos and K. Olive, |hep-ph/010233TD . 
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Figure 1: The (Tscaiar, as function of tan/?, for both positive and negative values 
of the parameter /x for given mo = 200 GeV, Ao = GeV. The sohd, dashed 
and dot-dashed hues correspond to M1/2 = 200, 400 and 600 GeV respectively. 
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Figure 2: Scattered plot of the dscaiar versus m^ from a sample of 5000 random 
points in the parameter space. Both cosmologically acceptable, Q^ h^ = 0.15 ± 
0.07 (marked as "In"), and unacceptable points (marked as "Out") are plotted. 
Low m^ values are excluded by chargino searches. The shaded region on the 
top is excluded by CDMS experiment M. 
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Figure 3: The same as in Figure |^, where only the cosmologically acceptable 
points are plotted. 
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Figure 4: The same as in Figure §, where the coannihilation corrections are 
taken into account. 
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Figure 5: Scattered plot of the (Tscaiar versus tan/? using the same random 
sample as in the Figure |. 
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